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Abstract 
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1 Introduction 



Self-sustained oscillations are found very often in nature. There are 
many examples in different branches of science such as in biology, chem- 
istry, mechanics, electronics, fluid dynamics, etc. [[Andronov et al., 1989|, 



Bai-Lin, 1990|| . Nonlinearities are required in order to have this kind of 
behaviour. The system reaches an oscillatory dynamics caracterized by 
a preferred period, wave form and amplitude, stable under slight pertur- 
bations. The oscillations are generated by an internal balance of ampli- 
fication and dissipation, even in the absence of external periodic forcing. 
(For instance, a nonlinear damping force which increases the amplitude 
for small velocities and decreases it for large velocities). This dynam- 
ical state can be modelled by the stable limit cycles found in specific 
nonlinear differential equations. 

Limit cycles are isolated closed trajectories in phase space. They are 
stable if the neighbouring solutions tend to them in an asymptotic sense 
or unstable if the neighbouring solutions unwind from them. Determi- 
nation of the number, amplitude and loci of limit cycles in a general 
nonlinear system is an unsolved problem that has attracted much atten- 
tion in this century ||Yan-Quian et al., 1986||. This constitues a part of 



the Hilbert's Sixteenth Problem ||Hilbert, 1902|| when we are restricted to 



two-dimensional autonomous systems of the form: 

x = P n (x,y), 
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y = Qn(x,y), 



(1) 



where P n and Q n are polynomials of degree n with real coefficients. Al- 
thoug it has been proved that the number of limit cycles in systems (f ) 
is finite | Ecallc ct al., 1987| , llyashcnko, f990 , the determination of the 



maximal number H n of limit cycles is still far away of being known. 

The van der Pol oscillator x + e(x 2 — l)x + x = 0, where x(t) = 
dx(t)/dt, is an example of system (1) that has been exhaustively studied. 
In this case, P 3 (x,y) = y and Q 3 (x,y) = —e(x 2 — l)y — x. It displays a 
limit cycle whose uniqueness and non-algebraicity has been shown for the 
whole range of the parameter e. Its behaviour runs from near-harmonic 
oscillations for e close to zero (e — > 0) to relaxation oscillations when e 
tends to infinity (e — > oo), making it a good model for many practical 
situations |[van der Pol, 1927| , [Lopez- Ruiz fc Pomeau, 1997|| . 



A generalization of the van der Pol oscillator is the Lienard equation, 

x + ef(x)x + x = 0, (2) 

with e a real parameter and f(x) any real function. When f(x) is a 
polynomial of degree iV = 2n + 1 or 2n this equation of the form (1) with 
Pn+i( x ,v) — V an d Qn+i(%, y) = ~ e f( x )y — x - It has been conjectured 
by Lins, Melo and Pugh (LMP-conjecture) that the maximum number 
of limit cycles allowed is just n [[Lins et al., 1977 1- It is true if iV = 2, or 



iV = 3 or if f(x) is even and iV = 4 ||Lins et al.,T977l |Rychkov, 19751 . 



Also, it is true in the strongly nonlinear regime (e — > oo) when f(x) 
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is an even polynomial [[Lopez fc Lopez- Ruiz, 200~0| . There are no gen 



eral results about the limit cycles when f(x) is a polynomial of degree 
greater than 5 neither, in general, when f(x) is an arbitrary real function 
[ Giacomini fc Neukirch, 1997 , piacomini fc Neukirch, 1998 ] . 



In the present paper, we are interested in the Lienard equation when 
f(x) is a continous even function, otherwise arbritary. We exploit the 
fact that the calculation of the number of limit cycles in the weak (e — > 0) 
and in the strongly (e — > oo) nonlinear regimes is possible for this kind 
of functions by means of simple algorithms. In fact, we find exactly the 
bifurcation curves of limit cycles in both regimes for several examples 
of viscous terms f(x). Section 2 is devoted to explain the strategies (or 
algorithms) used to calculate the amplitude and number of limit cycles 
in those extreme regimes, and in Section 3 we analyze some particular 
cases found in the literature. Last Section contains the conclusions. 



2 Limit Cycles in the Lienard Equation 

In order to study the limit cycles of equation (2) it is convenient to rewrite 
it in the coordinates (x,x) = (x,y) in the plane. We perform the change 
of variables x(t) = y(x) and x{t) = y(x)y'(x) (where y'(x) = dy/dx): 

yy' + ef(x)y + x = 0. (3) 

A limit cycle C\ = (x, y±{x)) of equation (3) has a positive branch 
y+{x) > and a negative branch y~(x) < 0. They cut the x-axis in 
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two points (— ai,0) and (02, 0) with 01,02 > 0. The oscillation x runs in 
the interval — ai < x < a<i- 

The origin (0, 0) is the only fixed point of equation (3). Then every limit 
cycle Ci solution of Eq. (3) encloses the origin. The result is a nested set 
of closed curves that defines the qualitative distribution of the integral 
curves in the plane (x,y). The stability of the limit cycles is alternated. 
For a given stable limit cycle, the two neighbouring limit cycles, the clos- 
est one in its interior and the closest one in its exterior, are unstable, and 
viceversa (Fig. 1). 

When f(x) is an even function, the symmetries of the equation (3) 
impose some properties over the shape of the limit cycles. Thus the 
inversion symmetry (x,y) <-> (—x,—y) implies y+(x) = — y~(— x) and 
ai = a 2 = a. Therefore, we can restrict ourselves to the positive branches 
of the limit cycles (x,y+(x)) with —a < x < a. The amplitude of os- 
cillation a identifies the limit cycle. The parameter inversion symmetry 
(e,x,y) <-> (— e, x, —y) implies that if C\ = (x,y±(x)) is a limit cycle for 
a given e, then Ci = (x, —y T (x)) is a limit cycle for — e. Moreover if 
C{ is stable (or unstable) then Ci is unstable (or stable, respectively). 
Therefore it is enough to consider the positive y-branch y+(x) of the 
limit cycles when e > for obtaining all the periodic solutions. (The 
limit cycles for a given — e < are obtained from a reflection over the 
x-axis of those limit cycles obtained for e > 0). 

Another property of a limit cycle can be derived from the fact that 
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the mechanical energy E = (x 2 + y 2 )/2 is conserved in a half oscillation: 



Thus, if equation (3) is integrated along the positive branch y+(x) of a 
limit cycle, between the maximal amplitudes of oscillation, we obtain: 



The solutions y+(x) of equation (3) and (4), vanishing in the extremes, 
constitute the finite set of limit cycles of equation (3). 

2.1 The Weakly Nonlinear Regime 

Lienard system (3) reduces to the simple harmonic oscillator when e = 0. 
All the circles y(x) = \Jr 2 — x 2 of radius r about the origin are solutions. 
This path-diagram is destroyed when e is slightly perturbed. Only the 
limit cycles survive as closed curves. They will have a slightly modified 
circular form. At order zero in e, we can suppose them as circles y+(x) = 
\J a 2 — x 2 with different amplitudes a's. Obviously, at this order, the 
condition (3) is verified, and condition (4) reads: 



Each solution ±a of the equation /5(a) =0 is the amplitude of a limit 
cycle of the Lienard system in the weak nonlinear regime. And viceversa, 
the amplitudes of all limit cycles of equation (3) are solutions of equation 
(5) in that regime. These results are exact for e = 0. In conclusion, 





(4) 
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equation (5) determines the amplitudes of the limit cycles of Lienard 
system defined by f(x) when e — > 0. 



The stability of a limit cycle in this regime is given, at the lowest 
order in e, by the sign of the integral: 



61() 



ef(x) 



e 

a 



dp (a) 
da 



J ao 



-ao y+\ x ) 

where a > is a solution of equation (5). The limit cycle is stable for 
a < and unstable for a > 0. 



As an example, we integrate equation (5) when f(x) is an even poly- 
nomial of degree 2n: 

f{x) = b + b 2 x 2 + + • • • + b 2n x 2n , 

where b , b 2 , 64, . . . , b 2n are real coefficients. Then, only the amplitudes a 
that satisfy the equation: 

Tra 2 " , (2A;)! 2k n 

8(a) = > b 2k -r-± — -rr-r;a = 

2 ^ 2k A k {k + 1)! k\ 

are allowed. The solution a = corresponds to the fixed point (0, 0) 

and the factor a 2 can be eliminated. Thus, the possible amplitudes a are 

the zeros of an even polynomial of degree In. There are no more than 

n different solutions a > and therefore, the maximun number of limit 

cycles in this case is n. We conclude that LMP-conjecture is true in the 

weak nonlinear regime. For instance, f(x) = x 2 — 1 in the van der Pol 

oscillator. Then P(a) = ixa 2 (a 2 — 4)/8 and the only existing limit cycle 

has the amplitude a ~ 2 when e — > 0. It is stable if e > and unstable 

if e < 0. 
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2.2 The Strongly Nonlinear Regime 



An algorithm that determines the number and amplitude of the limit 
cycles of Lienard systems in the strongly nonlinear regime has been pro- 
posed in reference [[Lopez fc Lopez- Ruiz, 2000f . A first approach to the 



shape of limit cycles when e —>■ +00 shows that the positive y-branch, 
y+(x) = ez?(x), of a stable limit cycle with amplitude a\ > is given by: 

s( \ _ / if - af < x < Si ( s 

Zi{X> \ -F(x) + F( Si ) if Si < x < af, { ' 

where F(x) = f$ f(t)dt and Sj < is called the gluing point of the 
two pieces z(x) = and z(x) = —F(x) + F(si). The unstable ones, 
y + (x) = ez"(a;), with amplitude a" > are given by: 

i{} \ if Ui<x<af, {) 

where Ui > is the gluing point of the two pieces in this case. 

In the remaining of this section we give a brief skecht of the algorithm 
(for a detailed discussion see reference ||L6pez fc Lopez- Ruiz, 2000| ). 



STABLE CYCLES : Consider the points s* < where F(x) has a positive 
local maximum and find the points a* defined by the rule: 

a* = min {x > s*, F(x) = F(s*)} . 

Geometrically a* represents the x-coordinate of the first crossing point 
between the straight z = F(s*) and the curve z = F(x) in the plane 
(x, z). If a* < \s*\ it is not possible to build the limit cycle and we can 
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eliminate this s* as a possible gluing point. If a* > \s*\ the point s* 
is a gluing point candidate. We rename and order all the pairs (s*,a*) 
verifying this last property as (sj, af) with s i+1 < Sj < and collect them 
into the set: 

A s = {(sj, af), F(si) local maximum, F(si) > 0,af > 

By construction af +1 > af. There are two different situations when two 
sucessive pairs, («i,af) and (sj+i,af +1 ), are considered: 

(a) — af +1 < Sj + i < — af < Sj. In this case it is possible to build a two- 
piecewise limit cycle with the pair (sj,af) as indicated in Eq. (6). This 
pair is picked out and renamed again as (sj, af). 

(b) — af +1 < —a* < Sj+i < Sj. Now the constuction of a limit cycle 
derived from the pair (sj, af) is not possible. This pair is rejected. 

If there is only one pair (si, af), we consider it satisfies (a). 

All the existing stable limit cycles can be found comparing the pairs 
i and i + 1 under rules (a)-(b) and iterating this process. All the pairs 
selected by condition (a) (and renamed as (sj,a|)) are collected into the 
set: 

A s = {( Si , a^)} = {{s h af) E A s , (s h a?) verifies (a)} . (8) 

The number, l s = card(A s ), of pairs (s,, af) is the number of stable limit 
cycles of the system (3). 



UNSTABLE CYCLES : The same process can be repeated for the 
unstable cycles by considering the points u* > 0, where F(u*) is a positive 
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local maximum, and their partners a* are defined by: 

a* = max {x < u*, F(x) = F(u*)} . 

The gluing point candidates u* must verify \a*\ > u* . After renaming 
and ordering the pairs (u*, a*) fulfilling this last condition as (uj, af) with 
u i+ i > Ui > 0, we collect them into the set: 

A u = {(Ui, af), F(Ui) local maximum, F(ui) > 0, |a"| > v-i} . 

A similar algorithm as indicated above can be applied in this case with 
the following modified rules: 

(a') Ui < —a" < Ui + \ < — a" +1 . In this case there exists an unstable 
two-piecewise limit cycle resulting from the pair (ui, af) and given in Eq. 
(7). This pair is picked out and renamed (u i: a"). 

(b') Ui < Ui + \ < —a" < —af +1 . The pair (Ui, af) does not produce a limit 
cycle and is rejected. 

If there is only one pair (ui, a") we consider it satisfies (a'). 

We iterate the process given by rules (a')-(b'). All the pairs selected 
by condition (a') are collected into the set: 

A u = {(u h af)} = {{u u a u i) e A u , (Ui,af) verifies (a')} . (9) 

The number, l u = card(A u ), of pairs (ui,af) is the number of unstable 
limit cycles of system (3). Obviously, l s — 1 < l u < l s + 1. 



It was claimed in [[Lopez fc Lopez-Ruiz, 2000| that the total number 
I of limit cycles of Eq. (3) in the strongly nonlinear regime is I = l s + l u , 
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where l s and l u are the number of stable and unstable limit cycles respec- 
tively. The amplitudes of these limit cycles are given by the numbers a* 
and af, respectively 



We remark also that each pair of zeros isc, of f(x) produces at most 
one limit cycle. If f(x) is an even polynomial of degree In there will 
be at most n limit cycles. Therefore, LMP-conjecture is also true in the 
strongly nonlinear regime. For instance, in the van der Pol oscillator, 
F(x) = —x + x 3 /3 has an unique positive local maximum at s — — 1. 
The amplitude a of the only existing limit cycle when e — > oo is given 
by the solution of the relation F(—l) = F(a), that is, a = 2. Its shape, 
y+(x) = ez(x), is (up to order e -2 ) given by: 




if - 2 < x < -1 

|(-x 3 + 3x + 2) if -1 < x < 2, 



3 Bifurcation Curves in some Examples 



We apply in this section the results of the former section to particular 
examples that have been studied by different authors in the literature. 

Example 1: f (x) = x 2n — 1, with n = 1, 2, 3, • • •. This case repre- 



sents a generalization of the van der Pol oscillator [[van der Pol, 1927 1 . It 

has only a limit cycle for the whole range of the parameter e. 

(i) e — > : The amplitude a n , n = 1,2,3, of this limit cycle in the 
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weakly nonlinear regime is the solution of the equation f3{a n ) = 0. We 
obtain: 

a n = 2 2? 



n\(n + 1)! 



\ (2n)! 

If n = 1 then ai = 2 (van der Pol system) and if n — > oo the result is 
aoo=l. That is, 1 < a n < 2 for all values of n. 

(ii) 6 — > oo : The calculation of the amplitude a n of this limit cycle in the 
strongly nonlinear regime requires to find the positive local maxima of 
F(x) = Jq f(t)dt. In this case, the only local maximum is at x = — 1 
with the value F(—l) = 2n/ (2n + 1). The amplitude a n is therefore given 
by: 

a n = F- 1 (-*?-) > 

In particular, ai = 2 and if n — » oo the amplitude is = 1. Also, in 
this case, 1 < a n < 2 for each value of n. 

(iii) < e < oo : Computer simulations show that the amplitude a n of 
the unique limit cycle is slighty perturbed in this regime. 

Example 2: f(x) = (x 2 — l)(x 2 — k), where k is a real parameter 



(Fig. 2). This sytem was studied by Lloyd in Ref. ||Lloyd, 1987|| . He 



showed that it has no periodic solutions if 1/5 < k < 5, and he suggested 
that there exist k* and k*, depending on e, such that there are two 
periodic solutions if < k < or k > k*, while there are none if 
k* < k < k*. Moreover, he finds that k* < (7 + v^45)/2 for some positive 
e. Here, the values of (fc*, k*) in the weakly, (k , k°), and in the strongly, 
(fcoo, k°°), nonlinear regimes are calculated. For intermediate values of e 
it is found numerically that koo < k* < ko and k° < k* < k°°. 
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(i) e — > : The amplitudes a±, solutions of the equation /9(a) = 0, of the 
limit cycles in this regime are the positive values of the expression: 



If k < there is only one limit cycle of amplitude a+. If k > the 
sign of the discriminant A = {k + l) 2 — 8A; determines the number of 
periodic solutions. The roots of A are: k = 3 — ~ 0.17157 and 
k° = 3 + V8 ~ 5.82842. Then, if A; < fc < A is negative and 
there is no limit cycle. If < k < k or k > k°, A is positive and the 
system has two periodic solutions of amplitudes a±. In k = a limit 
cycle of small amplitude bifurcates from the origin after an Andronov- 
Hopf bifurcation. In k = k Q and k = k° the two limit cycles appear or 
disappear by a saddle-node bifurcation. 

(ii) e — » oo : The limit cycles in this regime are determined by the the 
positive local maxima of F(x) = x 5 /5 — {k + l)x 3 /3 + kx. If k < there 
is only a positive local maximum at x + = —1 with the value F(—l) = 
(2 — 10/c)/15. The amplitude a + of this limit cycle is: 



If0<fc<l/5 there are two positive local maxima: one is at x + = — 1 



(3 — v / 5) 2 /4 ~ 0.14589. The amplitudes a± of these limit cycles are: 





and the other one is at x_ = \fk. The condition for the existence of two 
limit cycles is: F(— 1) > F(\/k). This is verified when < k < k^ = 
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If 1/5 < k < 1 there is only a positive local maximum at x — y/k for 
which -F _1 (x) has not an antiimage. Then there is not periodic solutions 
in that interval of k. The same behaviour is found when 1 < k < 5, but 
now the positive local maximum is located at x = 1. 
If k > 5 there are two positive local maxima localized at x + = —\fk and 
at x_ = 1. The condition for having periodic solutions reads: F(—y/k) > 
F(l). This condition holds when k > k°° = (3 + \/5) 2 /4 ~ 6.85410. The 
amplitudes a± of these limit cycles are: 

! / (2ife - 10)fct 

(3-|_ = F 




max <F 1 

I V 15 

If fcoo < k < k°° the system has no periodic solutions. 

(hi) 0<6<oo : Numerical computations of this system suggest that 

the curves fc*(e) and k*(e) behave as it is shown in Fig. 2. Thus we find 

that koo <k*<kQ and k° < k* < k°° for every real e. 

In summary, if k < k < k° there is no periodic solution, and, if < 

k < k^ or k > k°° the system has two limit cycles. Observe that k°° < 

(7 + \/45)/2, and the requirement k > (7 + y/45)/2 is not necessary for 



having two limit cycles for some e, as it was proposed in [Lloyd, 1987 
Some amplitudes a± are given in the following table: 
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k 


e = 


e = 1 


a 





V2 


1.41431 


a + 











a_ 


0.16 


1.19001 


1.19000 


a + 


0.16 


0.95072 


0.95172 


a_ 


5 








a + 


5 








a_ 


7 


3.29065 


3.29504 


a + 


7 


2.27410 


2.44764 


a_ 



Example 3: f(x) = 5x 4 — 3/ix 2 + 5, where /i and 5 are two real 
parameters (Fig. 3). This system is a generalization of Example 2. 
Several authors have studied the case 5=1: in Ref. [Kychkov, 1975 



Rychkov shows that this equation have two cycles when e > and /i > 
2.5; Alsholm has improved this result lowering the bound to fi > 2.3178 
(// > 2.3178 5*), and in ||0dani, 1996|| , Odani obtained a sharper result 
/i > v5- m ||Giacomini fc Neukirch, 1997] , Giacomini & Neukirch obtain 
a sequence of algebraic approximations, in the parameter plane (5, /x) for 
e = 1, to the bifurcation set B e=1 (S, fi) = 0, where the system undergoes 
a saddle-node bifurcation. Here, the bifurcation curves Bq(8, y) = and 
-Boo (5, A*) — in the weakly and in the strongly nonlinear regimes are 
calculated, respectively. We obtain B (S,fi) = 9/i 2 — 405 and -Boo (5, a0 = 
/i 2 — 55. For intermediate values of e, numerical simulations show that the 
curves B e (8,fi) = are localized between B (S,fj,) = and -Boo (5, /-0 = 
in such a way that if e<i > ei then B e2 (5, /i) = is between B ei (5, /i) = 
and -Boo (5, /-0 = 0. These results are in agreement with the earlier works 
cited above and let us a better understanding of the behaviour of this 
system for all the values of 5 and \i. 
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(i) e — > : The amplitudes a± of the limit cycles in this regime are the 
positive solutions of the equation /3(a) = 0, that is: 



If 5 < there is only one limit cycle of amplitude a + . If 5 > and 
fi < there are no real solutions for a±. If S > and /i > 0, the sign of 
the discriminant B (5,ii) = 9/x 2 — 405 determines the number of periodic 



has two limit cycles of amplitudes a±. At 5 = a limit cycle of small 
amplitude bifurcates from the origin after an Andronov-Hopf bifurcation. 
For the values (5, /i) where B (5,fi) = the two limit cycles appear or 
disappear by a saddle-node bifurcation. 

(ii) e — » oo : The positive local maxima of F(x) = x 5 — fix 3 + 5x must 
be found. If 5 < there is only a positive local maximum at x = 



-(^§r)*A/I6 where A = - 205. We have F(ar ) = -(3^ 2 + 

/lA — 20S)x /25. The amplitude a + of this limit cycle is: 
_ 1 \(3fi 2 + fiA- 205\/3fi + Ay] 




solutions. If 3/x < V^O 5 2 there is none and if 3/i > V^O 5 a the system 




If 5 > and 3/i < V20 F(x) has no local maxima and the system 
has no periodic solutions. If 8 > and 3/i > v^O 5^ there are two 
positive local maxima localized at x + = —(3/i + A)^/vT0 and x_ = 
(3/i — A)5/vT0. The condition for having two limit cycles is: F(x+) > 
F(x_). This follows if /i > v^o"^^. The amplitudes a± of the limit cycles 



are: 
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max 



F- 1 



-3/i 2 + /iA + 205 \ /3/i-A\2 



25 



V 10 



< 



If [i < V55? there are not periodic solutions. Then the bifurcation curve 
where the system undergoes the saddle-node bifurcation is defined by 
Boc(S,n) = - 55. 

(iii) <C e <C oo : Numerical simulations of this system suggest that the 
bifurcation curves B t (S,fj,) = are located between B (5,fi) = and 
.Boof^, /x) = 0, in such a way that if e 2 > e\ then B t2 (5, /i) = is located 
between B ei (5,fi) = and B^^S, /z) = 0. If yU*(e) is the solution of 



B € (5 , /i) = for a fixed <5 , then v/405 /3 < //(e) < V55o and /x*(e 2 ) > 
if e 2 > ex. 

As an example we give some values of the amplitudes a± of the periodic 
solutions for jj, — 1: 



5 


e = 


e= 1 


a 


-1 


v/2 


1.40990 


a + 


0.1 


1.02334 


1.02344 


a + 


0.1 


0.39087 


0.39090 


a_ 


0.3 








a + 


0.3 








a_ 



Example 4: f (x) = 7x 6 - 5(29 + b 2 )x 4 + 3(100 + 29b 2 )x 2 - 100b s 

where h is a real parameter (Fig. 4). Giacomini & Neukirch have investi- 
gated this system in Ref. ||Giacomini fc Neukirch, 1997|| . They find that 



the solutions of the equation F(x, b) = do not give the right qualitative 
amplitude-bifurcation diagram, where F(x, h) = x(x 2 — b 2 )(x 2 — 2 2 )(x 2 — 
5 2 ). In fact, the plot of the roots of F(x,b) = announce the presence 
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of a transtricital bifurcation near b = 2 and 6 = 5, and an Andronov- 
Hopf bifurcation at b = 0. Their method of algebraic approximations to 
the bifurcation curves shows that the supposed transcritical bifurcations 
are indeed saddle-node bifurcations. They conclude that the system can 
have one or three limit cycles. We confirm these results and stablish 
the correct amplitude-bifurcation curves a\(b), i = 1,2,3, for the three 
limit cycles, in the weakly, a^(b), and in the strongly, a°°(6), nonlinear 
regimes. The values of b for which the saddle-node bifurcations occur 
are calculated in those regimes. Numerical simulations show that for 
intermediate values of e the amplitude-bifurcation a\{b) curves and the 
values of b for which the saddle-node bifurcations occur are localized in 
the regions bounded by the curves a°(6) and af(b), % — 1, 2, 3. As in the 
former examples, the variation of parameter e does not introduce new 
qualitative information in the system, and only produces slight pertur- 
bations in the amplitude-bifurcation diagrams. 

(i) e — » : Mapple calculations allow us to solve the equation (3(a) = for 
the amplitudes a®(b), i = 1, 2, 3. The number of positive real solutions of 
that cubic equation in a 2 is determined by the the sign of the polynomial 
A(6) = -0.01784 b 8 + 1.15301 b 6 - 21.65794 b 4 + 132.559 b 2 - 189.45. 
If A < there are three limit cycles and if A > there is only one. If 
A = a saddle-node bifurcation arises in the system. The positive roots 
of A(6) are: 6° = 1.42636, b° 2 = 2.84148, b° 3 = 4.17545 and b° A = 6.08945. 
If < b < 6?, b° 2 < b < bl or b > b° 4 then A(6) < and there are three 
periodic solutions . If b\ < b < b® or 6° < b < 6° then A (6) > and there 
is only one. The amplitudes of the limit cycles for a given b are the cuts 
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of the line b — b with the curves a?(fc), a%{b) and a^ib) (see Fig. 4). 

(ii) e — » oo : Numerical calculation of show that the values of b for which 
the system undergoes a saddle- node bifurcation are: bf = 1.21, bf = 
3.49, bf = 3.95 and bf = 6.40. If < b < bf, bf < b < bf or b > bf 
the system has three limit cycles, and if bf < b < bf or bf < b < bf 
there is only one periodic solution (see Fig. 4). 

(iii) <^ e <^ oo : Numerical computations show that the amplitude- 
curves al(b), i = 1,2, 3, are localized in the narrow shaded region bounded 
by a®(b) and af{b) (Fig. 4a). Remark that the behaviour of the system 
suggests, once more, that the number of limit cycles do not increase when 
e increases (Fig. 4b). 



4 Conclusions 



Limit cycles are isolated periodic solutions of specific nonlinear differen- 
tial equations and can model self-sustained oscillations in nature. There 
are two difficult and connected problems in relation with limit cycles: the 
determination of bifurcation curves of these solutions in the parameter 
space and the determination of the maximal number of such solutions. 
In this work, we have exploited the possibility of calculating the bifur- 
cation curves of the limit cycles of Lienard equation x + ef(x)x + x = 
in the weakly (e — > 0) and in the strongly (e — > oo) nonlinear regimes 
when the viscous term f(x) is even. Firstly, these calculations allow us 
to improve the results existing in the literature for different examples in 
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these regimes. Secondly, the systems analyzed seem to follow the same 
pattern: the number of limit cycles does not increase when the non- 
linearity e increases. Moreover, the bifurcation curves for intermediate 
(0 < e <C oo) nonlinearity are always located between the bifurcation 
curves corresponding to the two extreme regimes. This means that al- 
though the variation of the nonlinearity e introduces an important mod- 
ification of the time scale and wave form of the oscillation, it perturbs 
slightly its amplitude. Only if two limit cycles have a very close ampli- 
tude for a given e, there exists the possibility of collapse of those limit 
cycles by a saddle- node bifurcation when |e| increases. If the system loses 
these two limit cycles it do not recover them for a stronger nonlinearity 
e. If the amplitudes of the limit cycles are separated enough for a given 
e, the number of periodic motions is conserved when e is varied. 
In particular, if we restrict ourselves to even-polynomial viscous forces, 
this behaviour suggests that Lins-Melo-Pugh conjecture on the number of 
limit cycles of Lienard systems is true. This is so because the conjecture 
is true in the weakly nonlinear regime and, according to the behaviour 
above explained, it should be true for any other regime. 
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Figure Captions 



Figure 1: A typical phase portrait of Eq. (3). The limit cycles of 
amplitudes a iy i = 1,2, enclose the origin and have the symmetry 
(x,y) <-> (— x, —y). Stable and unstable limit cycles alternate. 

Figure 2: Qualitative bifurcation diagram in the parameter plane 
(k, e) of system (2) for f(x) given in Example 2. In region I, where k < 0, 
there is only one periodic solution; in region II, where < k < fc*(e) or 
k > k*(e), there are two limit cycles, and in region 0> where fc*(e) < 
k < k*(e), there are none. On the line k — the system undergoes an 
Andronov-Hopf bifurcation and on the curves fc*(e) and k*(e) a saddle- 
node bifurcation. (Nomenclature in the text: fc*(0) = k , k*(0) = 
fc*(oo) = /coo and /c*(oo) = A;°°). 

Figure 3: The complete bifurcation diagram of system (2) for f(x) 
given in Example 3. The system has no periodic solutions in region 0> 
one limit cycle in region / and two limit cycles in region //. On the line 
5 = the system undergoes an Andronov-Hopf bifurcation, and on the 
curves B (5,fi) = 9fi 2 — 405 = and B^S,/!) = /i 2 — 55 = a saddle- 
node bifurcation arises for e = and e = oo, respectively. In region B 
are located all the bifurcation curves B e (5,/j,): the curve B €2 is located 
between B ei and -Boo if ei < €2. 



24 



Figure 4: Bifurcation curves of system (2) for f(x) given in Example 
4: (a) Amplitude bifurcation-diagram where a°(6) correspond to e = 
and a°°(b) to e = oo, % — 1,2,3. The amplitude-bifurcation a\{b) curves 
are located in the interior of the shaded region for every e. The number 
and amplitudes of limit cycles for b = b are the number and a-coordinates 
of the intersections between the curves a\(b) and the line b — b. (solid 
lines correspond to stable cycles and dashed lines to unstable ones for 
e > 0). 

(b) Qualitative bifurcation diagram in the parameter plane (b, e). In 
region I there is only a periodic solution and in region /// there are 
three limit cycles. 
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